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In these lectures I first give a motivation for investigations of in-medium properties of 
hadrons. I discuss the relevant symmetries of QCD and how they might affect the observed 
hadron properties. I then discuss at length the observable consequences of in-medium 
changes of hadronic properties in reactions with elementary probes, and in particular pho- 
tons, on nuclei. Here I put an emphasis on new experiments on changes of the a and u> 
mesons in medium. 



1. Introduction 

Studies of in-medium properties of hadrons are driven by a number of - partly connected - 
motivations. A first motivation for the study of hadronic in-medium properties is provided 
by our interest in understanding the structure of large dense systems, such as the interior 
of stars. This structure obviously depends on the composition of stellar matter and its 
interactions (for a recent review see [1]). 

The second motivation is based on the expectation that changes of hadronic properties in 
medium can be precursor phenomena to changes in the phase structure of nuclear matter. 
Here the transition to the chirally symmetric phase, that exhibits manifestly the symmetries 
of the underlying theory of strong interactions, i.e. QCD, is of particular interest. Present 
day's ultrarelativic heavy-ion collisions explore this phase boundary in the limit of high 
temperatures (T sa 170 MeV) and low densities. The other limit (low temperatures and 
high densities) is harder to reach, although the older AGS heavy-ion experiments and the 
planned CBM experiment at the new FAIR facility [2] may yield insight into this area. 
However, even in these experiments the temperatures reached are still sizeable (T w 120 
MeV). At very low temperatures the only feasible method seems to be the exploration 
of the hadronic structure inside ordinary nuclei, at the prize of a low density. Here the 
temperature is T = and the density at most equals the equilibrium density of nuclear 
matter, po- It is thus of great interest to explore if such low densities can already give 
precursor signals for chiral symmetry restoration. 

That hadrons can indeed change their properties and couplings in the nuclear medium 
has been well known to nuclear physicists since the days of the Delta-hole model that 
dealt with the changes of the properties of the pion and Delta-resonance inside nuclei [3]. 
Due to the predominant p-wave interaction of pions with nucleons one observes here a 
lowering of the pion branch with increasing pion-momentum and nucleon-density. This 
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effect can be seen in optical model analyses of pion scattering on nuclei, but the absorptive 
part of the 7r-nucleus interaction limits the sensitivity to small densities. More recently, 
experiments at the FSR at GSI have shown that also the mass of a pion at rest in the 
nuclear medium differs from its value in vacuum [4]. This is interesting since there are 
also recent experiments [5] that look for in-medium changes of the a meson, the chiral 
partner of the pion. Any comparison of scalar and pseudoscalar strength could thus give 
information about the degree of chiral symmetry restoration in nuclear matter. 

In addition, experiments for charged kaon production at GSI [6] have given some evi- 
dence for the theoretically predicted lowering of the K~ mass in medium and the (weaker) 
rising of the K + mass. State-of-the-art calculations of the in-medium properties of kaons 
have shown that the usual quasi-particle approximation for these particles is no longer 
justified inside nuclear matter where they acquire a broad spectral function [7,8]. 

At higher energies, at the CERN SPS and most recently at the Brookhaven RHIC, in- 
medium changes of vector mesons have found increased interest. This was due to two 
main reasons. First, these mesons are the first excitations of the QCD vacuum that are not 
protected by chiral symmetry, as the pions as Goldstone bosons are, and should therefore 
be more directly related to the chiral condensates. This expectation was triggered by the 
original work of Hatsuda and Lee [9] that, based on QCD sum rules, predicted a significant 
lowering of the vector meson masses with increasing density, the effect being as large as 
30% already at nuclear matter equilibrium density. As discussed in [10,11] this original 
prediction strongly depended on simplifying assumptions for the spectral function of the 
particles involved. When more realistic spectral shapes are used the QCD sum rule gives 
only certain restrictions on mass and width of the particles involved, but does not fix the 
latter; for that hadronic models are needed. In particular, for the p meson it turned out that 
the broadening is more dominant than a mass-shift [12]. 

The second reason for the interest in vector meson production stems from the fact that 
these mesons couple strongly to the photon so that electromagnetic signals can yield in- 
formation about properties of hadrons deeply embedded into nuclear matter. Indeed, the 
CERES experiment [13] has found a considerable excess of dileptons in an invariant mass 
range from w 300 MeV to w 700 MeV as compared to expectations based on the assump- 
tion of freely radiating mesons. 

This result has found an explanation in terms of a shift of the p meson spectral function 
down to lower masses, as expected from theory (see, e.g., [12,14-16]). However, the 
actual reason for the observed dilepton excess is far from clear. Both models that just 
shift the pole mass of the vector meson as well as those that also modify the spectral 
shape have successfully explained the data [17-19]; in addition, even a calculation that 
just used the free radiation rates with their - often quite large - experimental uncertainties 
was compatible with the observations [20]. There are also calculations that attribute the 
observed effect to radiation from a quark-gluon plasma [21]. 

The more recent experimental results on a change of the rho-meson properties in- 
medium obtained in an ultrarelativistic Au + Au collision by the STAR collaboration at 
RHIC [22] show a downward shift of the p-meson pole mass by about 70 MeV. Since also 
the p + p data obtained in the same experiment show a similar, though slightly less pro- 
nounced (-40 MeV), shift of the p-meson mass, phase-space distortions of the p-meson 
spectral shape may be at least partly responsible for the observed mass shift. In the heavy- 
ion experiment then a number of additional effects, mainly by dynamical interactions with 
surrounding matter, may contribute, but are hard to separate from the more mundane phase 
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space effects [23]. 

For both experiments quite different model calculations tend to explain the data, though 
often with some model assumptions. Their theoretical input is sufficiently different as to 
make the inverse conclusion that the data prove one or another of the proposed explanations 
impossible. In particular, if one is after some 'exotic' effect like a signal for the QGP, one 
has to make sure that one understands the 'classical' contributions to the observed signals 
very well. 

I have therefore already some years ago proposed to look for the theoretically predicted 
changes of vector meson properties inside the nuclear medium in reactions on normal 
nuclei with more microscopic probes [24,25]. Of course, the nuclear density felt by the 
vector mesons in such experiments lies much below the equilibrium density of nuclear 
matter, p , so that naively any density-dependent effects are expected to be much smaller 
than in heavy-ion reactions. 

On the other hand, there is a big advantage to these experiments: they proceed with the 
spectator matter being close to its equilibrium state. This is essential because all theoret- 
ical predictions of in-medium properties of hadrons are based on an equilibrium model 
in which the hadron (vector meson) under investigation is embedded in cold nuclear mat- 
ter in equilibrium and with infinite extension. However, a relativistic heavy-ion reaction 
proceeds - at least initially - far from equilibrium. Even if equilibrium is reached in a 
heavy-ion collision this state changes by cooling through expansion and particle emission 
and any observed signal is built up by integrating over the emissions from all these different 
stages of the reaction. 

Calculations of in-medium properties of hadrons thus necessarily rely on a number of 
simplifying assumptions, foremost being that of an infinite medium at rest in which the 
hadron under study is embedded. The properties so calculated are then, in a second step, 
being locally inserted into a time-dependent event simulation. In actual experiments these 
hadrons are observed through their decay products and these have to travel through the 
surrounding nuclear matter to the detectors. Except for the case of electromagnetic signals 
(photons, dileptons) this is connected with often sizeable final state interactions (FSI) that 
have to be treated as realistic as possible. For a long period the Glauber approximation 
which allows only for absorptive processes along a straight-line path has been the method 
of choice in theories of photonuclear reactions on nuclei. This may be sufficient if one is 
only interested in total yields of strongly absorbed particles. However, it is clearly insuffi- 
cient when one aims at, for example, reconstructing the spectral function of a hadron inside 
matter through its decay products. Rescattering and sidefeeding through coupled channel 
effects can affect the final result so that a realistic description of such effects is absolutely 
mandatory [26]. 

In this lecture note I first outline the fundamentals of chiral symmetry, its expected 
density dependence and its connection with actual observables. I then summarize results 
that we have obtained in studies of observable consequences of in-medium changes of 
hadronic spectral functions in reactions of elementary probes with nuclei. I demonstrate 
that the expected in-medium sensitivity in such reactions is as high as that in relativistic 
heavy-ion collisions and that in particular photonuclear reactions present an independent, 
cleaner testing ground for assumptions made in analyzing heavy-ion reactions. 
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2. Theory 



2.1 Chiral Symmetry 



A large part of the current interest in in-medium properties of hadrons comes from the 
hope to learn something about quarks in nuclei. More specifically, one hopes to see pre- 
cursors of a restoration of the original symmetries of the theory of strong interactions, i.e. 
QCD, which are spontaneously broken in our world. In this section of the lecture I briefly 
summarize the - for the present discussion - most relevant features of QCD. This material 
here draws heavily on my book [27] where more many details can be found. 
The QCD Lagrangian 



£qcd = -\G C ^ V + ]T [qjiyfDp - m f )q } ] 



(1) 



/ 



describes all the strong interactions of color-triplets q representing the quarks. Here the 
sum over c runs over the eight members of the color octet and is given by 

D,=d, + ig^Gl . (2) 

The G^j represent the eight vector fields describing the gauge bosons of the theory, called 
"gluons", because they transmit the binding forces and G^ in (1) denotes the correspond- 
ing field thensor. g is a dimensionless coupling constant (for further details see Chap. 15 
in [27]). 

If one assumes that the masses of the up- and down-quarks are nearly the same, then 
QCD possesses a global SU(2)y symmetry under the transformation 

q c (x)^q' c (x) =J2e- iEf -q c (x) . (3) 
/=i 

which acts onto the flavor degrees of freedom of the quark spinor given by 

u c (x) 

q c (x) = | d c (x) | , (4) 
s c {x) 

where u c , d c and s c represent the up, down and strange quark, respectively, with a color 
index c. The sum over the flavor index / in equation (3) runs only from 1 to 3, i.e. over the 
generators A-^ of the SU(2) subgroup of the full flavor-S'[/(3). As a consequence of the 
invariance of £qcd under the transformation (3), the flavor vector current 

Vf{x)=q cl ^q c (/ = 1,2,3) (5) 

is conserved on the quark level. The observed isospin symmetry of strong interactions 
and the corresponding vector current conservation (CVC) can thus be understood as a 
consequence of a symmetry of quark-quark interactions. 
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If, in addition, m u = rrid = 0, then there is also a symmetry of the QCD Lagrangian 
under 



q c (x)^q' c (x)=J2e-' n5£ i<lc(x) , (6) 
/=1 

so that the axial vector current 

is also conserved. In these considerations one uses the fact that the gluon fields are pure 
color fields and are thus not affected by the flavor generators A-' or by the Dirac algebra, 
e.g. by 75. 

It is clear that the chiral symmetry cannot be a manifest one on the hadronic level, 
since the hadrons are all massive. However, the Goldstone mechanism of spontaneous 
symmetry breaking can reconcile the nonzero masses of the hadrons with the observed 
(partial) conservation of the axial current. The axial SU(2) symmetry is then assumed to 
be realized in the Goldstone mode in which the Lagrangian, but not the state of the system, 
possesses the axial SU (2) symmetry. Therefore, it is obvious that in the real world also 
the chiral symmetry must be spontaneously broken. An order parameter for this symmetry 
breaking is given by the so-called chiral condensate 

(|««|> ± , (8) 

where the expectation value is taken in the vacuum state. The chiral condensate and thus 
the ground state can be chirally symmetric only if this condensate is zero as one can show 
by applying the transformation (6) (see also discussion below). 

A very simple estimate shows that the chiral condensate in the nuclear medium is in 
lowest order in density given by [16] 



(qq) med (p, T) « 1 - V hF ^ F l (qq) vac . (9) 




Here p s is the scalar density of the hadron h in the nuclear system and H h the so-called 
sigma-commutator that contains information on the chiral properties of h. The sum runs 
over all hadronic states. While (9) is nearly exact, its actual value is limited because neither 
the sigma-commutators of the higher lying hadrons nor their scalar densities are known. 
However, at very low temperatures close to the groundstate these are accessible for the 
nucleon. Here p s w p v ^ so that the condensate drops linearly with the nuclear (vector) 
density. Inserting numerical values for the physical constants in (9) gives 

(OT)mcd(p, 0) « - 0.3-£) (^) vac . (10) 

This drop of the chiral condensate with density can be understood in physical terms: 
with increasing density the hadrons with their chirally symmetric phase in their interior 
fill in more and more space in the vacuum with its spontaneously broken chiral symmetry. 
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Note that this is a pure volume effect; it is there already for a free, non-interacting hadron 
gas. One can obtain a very simple estimate for this effect by assuming that inside baryons 
the chirally symmetric phase prevails, as it does in bag models (for a review of such models 
see [27]). Assuming a baryon radius of w 0.8 fm one obtains Vjv ~ 2.1 fm 3 . This number 
has to be compared with the specific volume of one baryon at normal nuclear matter density 
(po = 0.15 fm~ 3 ): v = 0.15" 1 fm 3 w 6.7 fm 3 . This shows that at about v/V N w 3 times 
nuclear equilibrium density the baryons significantly overlap so that their chirally restored 
phase fills all space. 

The so-called Nambu-Jona-Lasinio (NJL) model has also often been used to determine 
the drop of the chiral condensate with density (and temperature) and to establish a link 
between the chiral condensate and physical masses. In this model the QCD Lagrangian 
is replaced by a Lagrangian that exhibits the chiral symmetry of QCD manifestly and in 
which the gluon fields are assumed to be integrated out. The interactions mediated by them 
are modeled by a point-interaction of the quarks. The NJL model in the two-flavor version 
is given by 



2 2 

C = iqj^q + G (qq) + (qi^rq) 



(11) 



While the authors originally formulated this model in terms of nucleon fields (quarks were 
not known then) we have given here the Lagrangian in terms of quark fields q. The La- 
grangian (11) represents a system of massless fermions interacting through a contact in- 
teraction (the term in square parentheses). The constant G is a coupling constant with the 
dimension of (mass) -2 ; it is assumed to be positive so that the self-interaction of the quark 
fields is attractive. Since the interaction is local it corresponds to a (5-function two-body 
potential between the quarks in a non-relativistic language. 

In Hartree-Fock mean field theory the single particle Hamiltonian can be obtained by 
a variation of the vacuum energy corresponding to the NJL Lagrangian with respect to q. 
The ensuing single particle Hamiltonian exhibits a mass term 

M q = ^ = -2G(qq)q . (12) 
dq 

The mass can be evaluated in a vacuum of free states 

M= -2G(0\qq\0) 

= 2G lv - = . (o) 

The sum over states in the Dirac sea diverges; it thus has to be regularized by introducing an 
appropriate cut-off A for the sum over momenta. Going over to an integral representation 
gives 

M = 4 G [ - = *P (14) 

4l<A y/p* + M 2 (2tt) 3 

It is obvious that (14) has M = as a solution, corresponding to a chirally symmetric 
vacuum (ground-) state. However, for large enough couplings also a solution with M > 
is possible. In this case the vacuum state no longer exhibits the chiral symmetry of the 
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Figure 1. Scaled quark mass m = M/A in the NJL model as a function of the scaled 
coupling constant g = A 2 /ty 2 G for vanishing Fermi momentum (from [27]) 



underlying theory, QCD; the symmetry is spontaneously broken. The value of the 'chiral 
condensate' (12) 



gives a measure for the amount of chiral symmetry breaking: it is zero if chiral symmetry is 
present in the vacuum and nonzero otherwise. Equation (14) has exactly the same structure 
as the gap equation in the BCS theory of superconductivity, used, for example, to obtain 
the pairing gap as a function of the pairing interaction strength. The mass M plays here 
the role of the the pairing gap A. 

It is a well-known property of the gap equation that the gap decreases - and ultimately 
disappears - due to the blocking effect when states above the Fermi-level become occu- 
pied, either through a thermal excitation of the system or through the presence of odd 
nucleons. Exactly the same phenomenon appears here: if also positive-energy quark states 
are occupied, either through a temperature in the system or through a non-vanishing baryon 
density, the integral on the rhs of ( 14) extends over these states as well; because of the prop- 
erties of the positive-energy spinors these states appear with an opposite sign and give a 
contribution to the mass 



In order to compensate for this negative contribution, with G fixed, the masses then have 
to become smaller until they have to vanish altogether, when the occupation of the positive 
energy states is increased more and more. Thus, with increasing density the quarks become 
massless again and chiral symmetry is restored. The same effect takes place when the 
temperature of the system is increased. In that case states above the Fermi momentum are 
partially occupied and this blocking leads to a decrease of the mass, just as in BCS theory. 




(15) 




(16) 
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2.2 QCD Sum Rules and Hadronic Models 



In the NJL model the dropping of the chiral condensate with density and/or temperature di- 
rectly causes a drop of the mass because both are linearly proportional. This is no longer the 
case in complex hadronic systems. How the drop of the scalar condensate there translates 
into observable hadron masses is not uniquely prescribed. The only rigorous connection 
is given by the QCD sum rules that relates an integral over the hadronic spectral function 
to a sum over combinations of quark- and gluon-condensates with powers of 1/Q 2 . The 
starting point for this connection is the electromagnetic current-current correlator 

U^(q) = i Jd 4 xe i<lx {Tj^x)j v (0)) , (17) 

where the current is expressed in terms of quark-field operators. This correlator can be 
decomposed [10] 

IW<Z) = q,q,R(q 2 ) - 9 ^U isotr (q 2 ) . (18) 
The QCD sum rule then reads 

oo 

R (Q)-Q2+^-—J ds {s + Q 2 )s d9) 



with Q 2 := -q 2 > and some subtraction constants £j. Here R OPE represents a Wilson's 
operator expansion (OPE) of the current-current correlator in terms of quark and gluon 
degrees of freedom in the space-like region. On the other hand, i? HAD (s) in (19) is the 
same object for time-like momenta, represented by a parametrization in terms of hadronic 
variables. Experimentally, it can be determined by hadron production in e+e~ reactions 
[28]. It is usually written in a form that exhibits the asymptotic, high momentum behavior 
of QCD in a manifest form: 



Imi? HAD (s) = e(s - s) Imi? RES (s) + 0(s - s ) ^- ( 1 + —) 



(20) 



where so denotes the threshold between the low energy region described by a spectral 
function for the lowest lying resonance, Imi? RES , and the high energy region described 
by a continuum calculated from perturbative QCD. The second term on the rhs of (20) 
represents the QCD perturbative result that survives when s — > oo. 

The dispersion integral in (19) connects time- and space-like momenta. Eq. (19) also 
connects the hadronic (rhs) with the quark (lhs) world. It allows - after the technical 
step of a Borel transformation - to determine parameters in a hadronic parametrization of 
i? HAD (s) by comparing the lhs of this equation with its rhs. 

The operator product expansion of i? OPE on the lhs involves a separation of momentum 
scales between the hard scale Q 2 and the soft, nonperturbative quark- and gluon conden- 
sates [10,1 1]. Using the measured, known vacuum spectral function for i? HAD allows one 
to obtain information about the condensates appearing on the lhs of (19). 

Turning this argument around one can model the density dependence of the quark con- 
densates on the lhs of (19) by using, for example, the density dependence obtained in the 
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NJL model. The QCDSR then gives information on the density dependence of the hadronic 
spectral function i? HAD on the rhs of (19). 

Since the spectral function appears under an integral the information obtained is, how- 
ever, not very specific. However, Leupold et al. have shown [10,11] that the QCDSR 
provides important constraints for the hadronic spectral functions in medium, but it does 
not fix them. Recently Kaempfer et al have turned this argument around by pointing out 
that measuring an in-medium spectral function of the u meson could help to determine the 
density dependence of the chiral condensate [29]. 

Thus models are needed for the hadronic interactions. The quantitatively reliable ones 
can at present be based only on 'classical' hadrons and their interactions. Indeed, in lowest 
order in the density the mass and width of an interacting hadron in nuclear matter at zero 
temperature and vector density p v are given by (for a meson, for example) 

m* 2 = to 2 - A-K^tf mN {q Q ,6 = 0) p v 
m*T* = mT° - 47r3/ roJV (<Zo, 9 = 0) p v . (21) 

Here f m N{qo,Q — 0) is the forward scattering amplitude for a meson with energy q on a 
nucleon. The width T° denotes the free decay width of the particle. For the imaginary part 
this is nothing other than the classical relation T* — T° = vap v for the collision width, 
where a is the total cross section. This can easily be seen by using the optical theorem. 

Actually evaluating mass and width from (21) requires knowledge of the scattering am- 
plitude which can only be obtained from very detailed analyses of experiments. The s- 
channel contributions to this scattering amplitude are determined by the properties of nu- 
cleon resonances and these are often not very well known yet. Here, resonance physics 
meets in-medium physics. 

Unfortunately it is not a-priori known up to which densities the low-density expansion 
(21) is useful. Post et al. [12] have recently investigated this question in a coupled-channel 
calculation of selfenergies. Their analysis comprises pions, ^-mesons and p-mesons as 
well as all baryon resonances with a sizeable coupling to any of these mesons. The authors 
of [12] find that already for densities less than O.bpo the linear scaling of the selfenergies 
inherent in (21) is badly violated for the p and the it mesons, whereas it is a reasonable ap- 
proximation for the r\ meson. Reasons for this deviation from linearity are Fermi-motion, 
Pauli-blocking, selfconsistency and short-range correlations. For different mesons differ- 
ent sources of the discrepancy prevail: for the p and r\ mesons the iterations act against 
the low-density theorem by inducing a strong broadening for the D 13 (1520) and a slightly 
repulsive mass shift for the 6*11 (1535) nucleon resonances to which the p and the 77 meson, 
respectively, couple. The investigation of in-medium properties of mesons, for example, 
thus involves at the same time the study of in-medium properties of nucleon resonances 
and is thus a coupled-channel problem. 

Note that such a picture, in which the selfenergies of hadrons are generated by inter- 
actions with the surrounding baryons, also encompasses the change of the chiral conden- 
sate in (9), obtained there for non-interacting hadrons. If the spectral function of a non- 
interacting hadron changes as a function of density, then in a classical hadronic theory, 
which works with fixed (free) hadron masses, this change will show up as an energy- 
dependent interaction and is thus contained in any empirical phenomenological cross sec- 
tion. 
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2.3 Coupled Channel Treatment of Incoherent Particle Production 

Very high nuclear densities (2 — 8po) an d temperatures T up to or even higher than w 
170 MeV can be reached with present day's accelerators in heavy-ion collisions. Thus, 
any density-dependent effect gets magnified in such collisions. However, the observed 
signal always represents a time-integral over various quite different stages of the collision 
- non-equilibrium and equilibrium, the latter at various densities and temperatures. The 
observables thus have to be modelled in a dynamic theory. In contrast, the theoretical input 
is always calculated under the simplifying assumption of a hadron in stationary nuclear 
matter in equilibrium and at fixed density. The results of such calculations are then used in 
dynamical simulations of various degrees of sophistication most of which invoke a quasi- 
stationary approximation. In order to avoid these intrinsic difficulties we have looked for 
possible effects in reactions that proceed closer to equilibrium, i.e. reactions of elementary 
probes such as protons, pions, and photons on nuclei. The densities probed in such reac- 
tions are always < po, with most of the nucleons actually being at about 0.5po- On the 
other hand, the target is stationary and the reaction proceeds much closer to (cold) equi- 
librium than in a relativistic heavy-ion collision. If any observable effects of in-medium 
changes of hadronic properties survive, even though the densities probed are always < po, 
then the study of hadronic in-medium properties in reactions with elementary probes on 
nuclei provides an essential baseline for in-medium effects in hot nuclear matter probed in 
ultra-relativistic heavy-ion collisions. 

With the aim of exploring this possibility we have over the last few years undertaken a 
number of calculations for proton- [30], pion- [31,32] and photon- [33] induced reactions. 
All of them have one feature in common: they treat the final state incoherently in a coupled 
channel transport calculation that allows for elastic and inelastic scattering of, particle pro- 
duction by and absorption of the produced hadrons. A new feature of these calculations is 
that hadrons with their correct spectral functions can actually be produced and transported 
consistently. This is quite an advantage over earlier treatments [34,35] in which the mesons 
were always produced and transported with their pole mass and their spectral function was 
later on folded in only for their decay. The method is summarized in the following section, 
more details can be found in [33]. 

We separate the photonuclear reaction into three steps. First, we determine the amount of 
shadowing for the incoming photon; this obviously depends on its momentum transfer Q 2 . 
Second, the primary particle is produced and third, the produced particles are propagated 
through the nuclear medium until they leave the nucleus. 

a Shadowing. Photonuclear reactions show shadowing in the entrance channel, for 
real photons from an energy of about 1 GeV on upwards [36]. This shadowing is due to 
a coherent superposition of bare photon and vector meson components in the incoming 
photon and is handled here by means of a Glauber multiple scattering model [26]. In 
this way we obtain for each value of virtuality Q 2 and energy v of the photon a spatial 
distribution for the probability that the incoming photon reaches a given point; for details 
see [26,37,38]. 

b Initial Production. The initial particle production is handled differently depending 
on the invariant mass W — y/s of the excited state of the nucleon. If W < 2 GeV, 
we invoke a nucleon resonance model that has been adjusted to nuclear data on resonance- 
driven particle production [33] . If W > 2 GeV the particle yield is calculated with standard 
codes developed for high energy nuclear reactions, i.e. FRITIOF or PYTHIA; details are 
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given in [39]. We have made efforts to ensure a smooth transition of cross sections in the 
transition from resonance physics to DIS. 

c Final State Interactions. The final state is described by a semiclassical coupled chan- 
nel transport model that had originally been developed for the description of heavy-ion 
collisions and has since then been applied to various more elementary reactions on nuclei 
with protons, pions and photons in the entrance channel. 

In this method the spectral phase space distributions of all particles involved are prop- 
agated in time, from the initial first contact of the photon with the nucleus all the way to 
the final hadrons leaving the nuclear volume on their way to the detector. The spectral 
phase space distributions Fh(x,p, fj,, t) give at each moment of time and for each particle 
class h the probability to find a particle of that class with a (possibly off-shell) mass p, and 
momentum p at position x. Its time-development is determined by the BUU equation 

(1 + ^£ - ^Jk)FH = G h A h L h F h . (22) 
at op or or op 

Here Hh gives the energy of the hadron h that is being transported; it contains the mass, the 
selfenergy (mean field) of the particle and a term that drives an off-shell particle back to its 
mass shell. The terms on the lhs of (22) are the so-called drift terms since they describe the 
independent transport of each hadron class h. The terms on the rhs of (22) are the collision 
terms; they describe both elastic and inelastic collisions between the hadrons. Here the 
term inelastic collisions includes those collisions that either lead to particle production or 
particle absorption. The former is described by the gain term GhAh on the rhs in (22), the 
latter process (absorption) by the loss term L h F h . Note that the gain term is proportional 
to the spectral function A of the particle being produced, thus allowing for production of 
off-shell particles. On the contrary, the loss term is proportional to the spectral phase space 
distribution itself: the more particles there are the more can be absorbed. The terms Gh 
and Lh on the rhs give the actual strength of the gain and loss terms, respectively. They 
have the form of Born-approximation collision integrals and take the Pauli-principle into 
account. The free collision rates themselves are taken from experiment or are calculated 
[33]. 

Eq. (22) contains a selfconsistency problem. The collision rates embedded in G and 
L determine the collisional broadening of the particles involved and thus their spectral 
function A. The widths of the particles, resonances or mesons, thus evolve in time away 
from their vacuum values. In addition, broad particles can be produced off their peak mass 
and then propagated. The extra 'potential' in H already mentioned ensures that all particles 
are being driven back to their mass-shell when they leave the nucleus. The actual method 
used is described in [33]. It is based on an analysis of the Kadanoff-Baym equation that has 
led to practical schemes for the propagation of off-shell particles [40,41]. The possibility 
to transport off-shell particles represents a major breakthrough in this field. For further 
details of the model see Ref. [33] and [39] and references therein. 
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3. Particle Production on Nuclei - Observables 



4. i] Production 

We first look at the prospects of using reactions with hadronic final states and discuss the 
photo-production of ?7-mesons on nuclei as a first example. These mesons are unique in 
that they are sensitive to one dominating resonance the Su (1535) so that one may hope to 
learn something about the properties of this resonance inside nuclei. Experiments for this 
reaction were performed both by the TAPS collaboration [42,43] and at KEK [44]. 

Estimates of the collisional broadening of the 5n(1535) resonance have given a colli- 
sion width of about 35 MeV at p [45]. The more recent, and more refined, selfconsistent 
calculations of [12] give a very similar value for this resonance. In addition, a dispersive 
calculation of the real part of the selfenergy for the resonance at rest gives only an insignif- 
icant shift of the resonance position. Thus any momentum-dependence of the selfenergy as 
observed in photon-nucleus data can directly be attributed to binding energy effects [46]. 
The results obtained in [46] indicate that the momentum dependence of the TV* (1535) 
potential has to be very similar to that of the nucleons. 

The particular advantage of our coupled-channels approach can be seen in Fig. 2 which 
shows results both for photon- and electroproduction of r]'s on nuclei; for the latter process 
so far no data are available. The calculations give the interesting result that for photo- 
production a secondary reaction channel becomes important at high energies and virtuali- 
ties. In this channel first a pion is produced which travels through the nucleus and through 
ttN — ► iV*(1535) — > Nt] produces the 77 that is finally observed in experiment. This 
channel becomes even dominant for high Q 2 electroproduction. In both cases the reason 
for the observed growing of the importance of secondary production channels lies in the 
higher momentum transferred to the initial pion [47]. 

4.1 2n Production 

If chiral symmetry is restored, the masses of the scalar isoscalar a meson and that of the 
scalar isovector pion should become degenerate. This implies that the spectral function 
of the a should become softer and narrower with its strength moving down to the lit 
threshold. This leads to a threshold enhancement in the nn invariant mass spectrum due to 
suppression of the phase-space for the a — > tttt decay. 

A first measurement of the two pion invariant mass spectrum has been obtained by 
the CHAOS collaboration in pion induced reactions on nuclei [48]. The authors of [48] 
claimed to indeed have seen an accumulation of spectral strength near the 2tt threshold for 
heavy target nuclei in the ir + -K~ mass distribution. According to the arguments presented 
in the introduction, photon induced reactions in nuclei are much better suited to investigate 
double pion pion production at finite baryon densities. A recent experiment with the TAPS 
spectrometer at the tagged-photon facility MAMI-B in Mainz indeed shows an even more 
pronounced accumulation of spectral strength of the two pion mass spectrum for low in- 
variant masses with increasing target mass corresponding to increasing average densities 
probed [5]. This accumulation has been observed for the 7r°7r° but not for the ^ir final 
state. 
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Figure 2. Eta-electroproduction cross sections on 40 C a for three values of Q 2 given 
in the figures. The uppermost curve in all three figures gives the total observed yield, 
the next curve from top gives the contribution of the direct r\ production channel and 
the lowest curve shows the contribution of the secondary process nN — > r/N. For the 
highest value of Q 2 = 3.6 GeV 2 the primary and secondary production channels nearly 
coincide (from [47]). 



These results have been explained by a model developed by Roca et al. [49]. In this 
model the a meson is generated dynamically as a resonance in the tttt scattering amplitude. 
By dressing the pion propagators in the medium by particle-hole loops, they found the 
observed downward shift of the tttt mass spectrum to be consistent with a dropping of the 
er-pole in the tttt scattering amplitude, i. e. a lowering of the a meson mass. 

In analyzing the TAPS data it is absolutely essential to simulate the final state interac- 
tions of the outgoing pions, correlated or not, in an as realistic way as possible. We have, 
therefore, taken the conservative approach of analyzing this reaction without any changes 
of the (7 spectral function in the outgoing channel. The result of this study [50] is shown in 
Fig. 3. 

Fig. 3 shows on its left side that the observed downward shift of the 2tt° mass spectrum 
can be well reproduced by final state interactions on the independent pions without any in- 
medium modification of the tttt interaction. This shift can be attributed to a slowing down 
of the pions due to quasi elastic collisions with the nucleons in the surrounding nuclear 
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Figure 3. Two pion invariant mass distributions for the (left) and 7T°7r (right) 
photoproduction on 12 C and 208 Pb (from [50]). The solid lines represent our results 
whereas the dashed lines labelled "Valencia" depict the results of ref. [49]. 



medium. Also shown is the result of the calculations of Roca et al. [49]. Both calculations 
obviously agree with each other so that the final observables are fairly independent of any 
7T7T correlations. 

On the right side of Fig. 3 the results for the semi-charged ir ^ channel are shown. 
Here our theoretical result (solid line) overestimates the data by up to a factor of 3. The re- 
sult obtained in Ref. [49] also lies too high, but by a lesser amount; the difference between 
the two calculations has been explained in [50] by the neglect of charge-transfer channels 
in the calculations of Roca et al. [49]. The observed discrepancy with experiment for this 
channel is astounding since the method used normally describes data within a much nar- 
rower error band. Thus understanding this discrepancy is absolutely essential before a shift 
or non-shift in the mass distribution for this channel can be ascertained. A new analysis of 
the data may be helpful in this regard [51]. 

Even before this problem has been resolved it is clear now that pion rescattering, in- 
cluding charge transfer, plays a major role for the 2tt mass distributions. It is, therefore, 
absolutely mandatory that any theoretical evaluation of this reaction takes this FSI effect 
into account. Because the reaction is incoherent quantum mechanical approaches to this 
problem are not feasible. The semiclassical transport approaches, on the other hand, suffer 
from intrinsic problems at very low pion energies (in the present experiment « 10 MeV 
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kinetic energy). Here, first, the de Broglie wavelength of the pion is large so that semi- 
classical approaches have to be critically examined. Second, for these low energies the 
low-energy real pion potential becomes important and has to be included together with the 
absorptive part. Careful analyses and comparisons with optical model approaches have 
shown [52] that both effects introduce some systematic error into the calculations that have 
to be taken into account in drawing any conclusions from this experiment. 

4.2 lo Production 

Many of the early studies of hadronic properties in medium concentrated on the p-meson 
[14,53], partly because of its possible significance for an interpretation of the CERES ex- 
periment. It is clear by now, however, that the dominant effect on the in-medium properties 
of the p-meson is collisional broadening that overshadows any possible mass shifts [12] and 
is thus experimentally hard to observe. The emphasis has, therefore, shifted to the u> me- 
son. An experiment measuring the A(-f, lo — ► tt q j')X reaction is presently being analyzed 
by the TAPS/Crystal Barrel collaborations at ELS A [54]. The varying theoretical predic- 
tions for the to mass (640-765 MeV) [53] and width (up to 50 MeV) [31,55] in nuclear 
matter at rest encourage the use of such an exclusive probe to learn about the lo spectral 
distribution in nuclei. 




Figure 4. Mass differential cross section for photoproduction off Nb. Shown 
are results both with and without a mass-shift as explained in [57]. The quantity p gives 
the escape probability for one of the four photons in the 2tv° channel. The two solid 
curves give results of calculations with p = 26 % with and without mass-shift, the two 
dashed curves give the same for p = 5% (from [56] 

Simulations have been performed at 1.2 GeV and 2.5 GeV photon energy, which cover 
the accessible energies of the TAPS/Crystal Barrel experiment. After reducing the combi- 
natorial and rescattering background by applying kinematic cuts on the outgoing particles, 
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we have obtained rather clear observable signals for an assumed dropping of the u> mass 
inside nuclei [57]. Therefore, in this case it should be possible to disentangle the collisional 
broadening from a dropping mass. 

One of the problems in the experimental analysis is the proper subtraction of the back- 
ground observed [58]; this has so far been done in a heuristic way [59]. Our calculations 
represent complete event simulations. It is, therefore, possible to calculate these back- 
ground contributions and to take experimental acceptance effects into account. An exam- 
ple is shown in Fig. 4 which shows the effects of a possible misidentification of the ui 
meson. This misidentification can come about through the 2ir° —> 47 channel if one of 
the four photons escapes detection and the remaining three photons are identified as stem- 
ming from the 7r°7 — > 37 decay channel of the w-meson. The calculations show that the 
misidentification does not affect the low-mass side of the omega spectral function. 

Fig. 4 shows a good agreement between the data of the TAPS/CB@ELSA collaboration 
[59] for a photon escape probability of 5 % and a mass shift = m° — 0.18 p/po- In 
[57] we have also discussed the momentum-dependence of the w-selfenergy in medium 
and have pointed out that this could be accessible through measurements which gate on 
different three-momenta of the ui decay products. The data presented by D. Trnka at this 
workshop [58] confirm this expectation. 

4.3 Dilepton Production 

Dileptons, i.e. electron-positron pairs, in the outgoing channel are an ideal probe for in- 
medium properties of hadrons since they experience no strong final state interaction. A 
first experiment to look for these dileptons in heavy-ion reactions was the DLS experiment 
at the BEVALAC in Berkeley [60]. Later on, and in a higher energy regime, the CERES 
experiment has received a lot of attention for its observation of an excess of dileptons with 
invariant masses below those of the lightest vector mesons [13]. Explanations of this excess 
have focused on a change of in-medium properties of these vector mesons in dense nuclear 
matter (see e.g. [17,18]). The radiating sources can be nicely seen in Fig. 5 that shows the 
dilepton spectrum obtained in a low-energy run at 40 AGeV together with the elementary 
sources of dilepton radiation. 

The figure exhibits clearly the rather strong contributions of the vector mesons - both 
direct and through their Dalitz decay - at invariant masses above about 500 MeV. The 
strong amplification of the dilepton rate at small invariant masses M caused by the photon 
propagator, which contributes <~ 1 /M 4 to the cross section, leads to a strong sensitivity 
to changes of the spectral function at small masses. Therefore, the excess observed in 
the CERES experiment can be explained by such changes as has been shown by various 
authors (see e.g. [34] for a review of such calculations). 

In view of the uncertainties in interpreting these results discussed earlier we have stud- 
ied the dilepton photo-production in reactions on nuclear targets. Looking for in-medium 
changes in such a reaction is not a priori hopeless: Even in relativistic heavy-ion reactions 
only about 1/2 of all dileptons come from densities larger than 2p [34]. In these reactions 
the pion-density gets quite large in the late stages of the collision. Correspondingly many 
p mesons are formed (through n + tt — > p) late in the collision, where the baryonic matter 
expands and its density becomes low again. 

In [33] we have analyzed the photoproduction of dileptons on nuclei in great detail. 
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Figure 5. Invariant dilepton mass spectrum obtained with the CERES experiment in 
Pb + Au collisions at 40 AGeV (from [13]). The thin curves give the contributions of 
individual hadronic sources to the total dilepton yield, the fat solid (modified spectral 
function) and the dash-dotted (dropping mass only) curves give the results of calcula- 
tions [19] employing an in-medium modified spectral function of the vector mesons. 



After removing the Bethe-Heitler contribution the dilepton mass spectrum in a 2 GeV 
photon-induced reaction looks very similar to that obtained in an ultrarelativistic heavy- 
ion collision (Fig. 5). The radiation sources are all the same in both otherwise quite dif- 
ferent reactions. The photon-induced reaction can thus be used as a baseline experiment 
that allows one to check crucial input into the simulations of more complicated heavy-ion 
collision. 

A typical result of such a calculation for the dilepton yield - after removing the Bethe- 
Heitler component - is given in Fig. 6. The lower part of Fig. 6 shows that we can expect 
observable effect of possible in-medium changes of the vector meson spectral functions 
in medium on the low-mass side of the u> peak. In [33] we have shown that these effects 
can be drastically enhanced if proper kinematic cuts are introduced that tend to enhance 
the in-medium decay of the vector mesons. There it was shown that in the heavy nucleus 
Pb the w-peak completely disappears from the spectrum if in-medium changes of width 
and mass are taken into account. The sensitivity of such reactions is thus as large as that 
observed in ultrarelativistic heavy-ion reactions. 

An experimental verification of this prediction would be a major step forward in our 
understanding of in-medium changes. The ongoing g7 experiment at JLAB is presently 
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compare with Fig. 5 (from [56]). 



analyzing such data [61]. This experiment can also yield important information on the 
time-like electromagnetic formfactor of the proton and its resonances [62] on which little 
or nothing is known. 

5. Conclusions 

In this lecture note I have first outlined the theoretical motivation for studies of in-medium 
properties of hadrons and their relation to QCD. I have then shown that photonuclear re- 
actions on nuclei can give observable consequences of in-medium changes of hadrons that 
are as big as those expected in heavy-ion collisions which reach much higher energies, 
but proceed farther away from equilibrium. Information from photonuclear reactions is 
important and relevant for an understanding of high density - high temperature phenom- 
ena in ultrarelativistic heavy-ion collisions. Special emphasis was put in these lectures not 
so much on the theoretical calculations of hadronic in-medium properties under simpli- 
fied conditions, but more on the final, observable effects of any such properties. I have 
discussed that for reliable predictions of observables one has to take the final state inter- 
actions with all their complications in a coupled channel calculation into account; simple 
Glauber-type descriptions are not sufficient. 

As an example that is free from complications by FSI I have shown that in photonuclear 
reactions in the 1 - 2 GeV range the expected sensitivity of dilepton spectra to changes 
of the p- and uj meson properties in medium is as large as that in ultrarelativistic heavy- 
ion collisions and that exactly the same sources contribute to the dilepton yield in both 
experiments. While the dilepton decay channel is free from hadronic final state interactions 
this is not so when the signal has to be reconstructed from hadrons present in the final state. 
While the ui photoproduction, identified by the semi-hadronic 7r°7 decay channel, seems 
to exhibit a rather clean in-medium signal, the double-pion experiments set up to look for 
an in-medium shift of scalar strength in nuclei are dominated by final state interactions of 
the produced pions that have to be properly taken into account in any reliable analysis. 
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